Propulsion with a Rotating Elastic Nanorod. 
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The dynamics of a rotating elastic filament is investigated using Stokesian simulations. The 
filament, straight and tilted with respect to its rotation axis for small driving torques, undergoes 
at a critical torque a strongly discontinuous shape bifurcation to a helical state. It induces a 
substantial forward propulsion whatever the sense of rotation: a nanomechanical force-rectification 
device is established. 
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Propulsion of a micrometer sized object through a vis- 
cous fluid, where friction dominates and motion is invari- 
ably overdamped, calls for design strategies very differ- 
ent from the ordinary macroscopic world, where inertia 
plays an important role . Biology has come up with a 
large number of concepts, ranging from active polymer- 
ization of gel networks yl, molecular motors moving on 
railway tracks formed by protein filaments , to rotating 
or beating propellerlike appendices Q- Those propellers 
involve stiff polymers that are moved by molecular mo- 
tors and rely on two different basic designs: i) Flagella of 
bacteria (prokaryotic cells) are helical stiff polymers set 
in motion at their base by a rotary motor. Hydrodynamic 
friction converts the rotational motion of the helix into 
thrust along the helix axis 1,4]. ii) Cilia and flagella of 
sperms (eukaryotic organisms) are rodlike polymers that 
are anchored to a surface and beat back and forth driven 
by internal motors. They are used by small cells to swim 
but also by internal organs to pump liquids 

For a number of biomedical applications, e.g. for di- 
rected motion of artificial viruses through cells or nano- 
devices through the bloodstream, as well as for mix- 
ing strategies in nanochips, it is desirable to develop 
similar synthetic propulsion mechanisms or to incorpo- 
rate biological single-molecule motors into synthetic en- 
vironments. Recent discoveries opened the route to the 
synthetic manufacture of rotary single-molecule motors 
driven by chemical ^ or optical HiQ energy. An ATPase 
can also be fixed to different substrates and used to rotate 
metallic or organic nanorods foj . This raises the ques- 
tion about the minimal design necessary to convert the 
rotational power of such nanoengines into directed thrust 
in a viscous environment. In this Letter, we show that an 
elastic rod that rotates around a point at a constrained 
azimuthal angle gives rise to a substantial forward thrust 
regardless of the sense of its rotation. It thus acts as a 
rectification device and produces net thrust even when it 
is stochastically rotated back and forth. It makes usage of 
a helical polymer unnecessary, allowing a selection from 
the much wider class of straight stiff polymers. A few 
theoretical works have elucidated the coupling between 
hydrodynamics in a viscous medium and elasticity of soft 
polymers [Hill El- It was shown that finite stiffness 
of beating straight filaments breaks time-reversal sym- 



metry and enables propulsion [Tol Il2j |. Conversely, re- 
versing the rotation sense of bacterial flagella gives rise 
to pronounced shape polymorphic transformations 
and shows that elasto-hydrodynamics are important in 
biology. Likewise, hydrodynamic interactions between 
rotating flagella are involved in bundling and are impor- 
tant for switching between tumbling and running modes 
of E. coli bacteria 0, . 

We study an elastic rod whose base rotates on a cone 
driven by an external torque applied to its anchoring 
point. We use Stokes-simulation techniques where hy- 
drodynamics are treated on the Rotne-Prager Green's 
function level including the full coupling between ther- 
mal, elastic and hydrodynamic forces. We obtain the 
nonlinear relation between driving torque and angular 
frequency in the presence of an external force load and 
estimate the power conversion efficiency. The filament 
exhibits a strongly discontinuous shape bifurcation at a 
finite torque value at which the rod jumps closer to the 
rotation axis. The elastic polymer is modeled as a chain 
of M -I- 1 connected spheres and the time evolution of 
the sphere positions ri{t) is governed by the position- 
Langevin equation 
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where the potential U{rk) = Y^fLo^ii^i^n+i - 2a)^ + 
^(1 — cos^i)] is the discrete version of the extensible 
wormlike chain model, a the sphere radius, r.y = — r^, 
and 9i is the angle between neighboring bonds of sphere i. 
The parameters 7 and e are the stretching and bending 
moduH, which are related by e/7 — (isotropic rod). 
Hydrodynamic interactions between two spheres i and j 
are implemented via the Rotne-Prager mobility tensor 
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which is an approximation valid for large sphere separa- 
tions [1^ . 1 denotes the unit 3x3 matrix and 77 the 
viscosity of the aqueous solvent. For the self-mobility 
we choose the Stokes mobility of a sphere fin = /io = 
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FIG. 1; a) Sketch in 3D and projected views of the rotating filament before (A, N = 600) and after (B, A'^ — 800) the bifurcation 
(M — 30 beads, £p/L = 10'^). b) Angular velocity uj across the shape transition as a function of time {N — 1.33 x 10'', 
£p/I/ — 2 X W^). c) Angular velocity as function of external torque, iV, for various £p/L: 3333 (o), 6667 (□), 1.33 x lO** (o) 
and 2 x 10'*(o) [inset: 333 (o), 1000 (□) and 2500 (o)]. d) Hysteresis cycle for £p/L = 1667 for increasing (o) and decreasing 
torque (□) at constant rate SN/{St/A) = 1.25 x 10~®. e) Persistence length vs. critical torque for increasing (o) and decreasing 
torque (□) both obeying a linear law according to Eq. Q (All data for stalled case, 14 = 0). 



l/(67r?7a). To mimic a rotary motor, we apply an ex- 
ternal force F™* on monomer 2, which is related to the 
torque N by F™* = N x ri2/rj2- The Langevin random 
displacement ^(i) models the action of the heat bath and 
obeys the fluctuation-dissipation relation {ii{t)^j{t')') = 
2kBTHi^S(t — t') numerically implemented by a Cholesky 
decomposition [l^. The persistence length is given by 
£p = e/{kBT). Twist and torsional degrees are omitted 
since for most synthetic polymers, free rotation around 
the polymer backbone is possible. Two different force 
ensembles are investigated: the stalled case, where the 
first two monomers forming the base are held fixed in 
space by applying virtual forces which exactly cancel all 
other elastic and hydrodynamic forces. In the moving 
case, we let the polymer move along the z axis and ap- 
ply virtual forces on the first two monomers only later- 
ally such that the polymer base moves along a vertical 
line. Hydrodynamic boundary effects are considered fur- 
ther below. A finite tilt angle at the base is imposed 
by a spontaneous curvature term in the elastic energy 
(e/2a)[l — cos(0i — 01)]; in this paper we show data for 
spontaneous curvature 6** = 45°. For the numerical it- 
erations, we discretize Eq. ^ with time step A. By 
rescaling time, space and energies, we arrive at rescaled 
parameters A = AfcsT^o/a^, ^p/L = s/{2aMkBT) 
and N = N/ksT. For sufRcient numerical accuracy we 
choose time steps in the range A = 10~'^-10^^. Output 
values are calculated every lO'^-lO'* steps, total simula- 
tion times are 10^ — 10^ steps. 



When a torque is applied to the filament base, it ro- 
tates and after some time exhibits a stationary shape. 
Since the friction is larger at the free end, the filament 
bends into a curved structure. The stationary angular ve- 
locity Co = bjo^ / {ksT Ho) (in degrees) is plotted in Fig.^ 
as a function of the applied torque N showing a nonlin- 
ear increase, which is caused by the finite bending rigid- 
ity. Indeed, for an infinitely stiff polymer, this relation 
is strictly linear. For a critical torque, A^c, a shape bi- 
furcation occurs and the angular frequency jumps dra- 
matically. Figure ^ shows the stationary shape of the 
polymer before (A) and after the bifurcation (B). The 
high-torque state is characterized by a smaller distance 
from the rotation axis. Figure^ shows the actual time- 
dependent behavior of the rotational speed as the rod 
crosses the transition and settles in the stationary state. 
Pronounced hysteresis (which becomes slightly weaker 
with decreasing torque sweep rate) is observed when N is 
varied across the critical region, as illustrated in Fig. 
Figure ^5 shows the dependence of the critical torque 
on the persistence length of the rod. Since the torque 
change rate is finite, increasing (o) and decreasing torque 
(□) gives rise to slightly different bifurcation values. For 
very small persistence lengths, conformational fluctua- 
tions suppress the transition. The shape bifurcation can 
be understood by simply balancing elastic and driving 
torques: the bending torque due to the filament deforma- 
tion, projected along the vertical axis, reads ~ esmOi/R 
where R is the bending radius of the filament, which has 
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FIG. 2: a) Propulsion velocity parallel to the rotation axis vs. NL/lp in the z-moving case for £p/L = 2500 (x), 5000 (□) 
and 2 x 10* (o). The inset shows the variations of the ratio of velocities after and before the transition vs. the ratio of angular 
velocities (o, £p/L — 1667). b) Vz vs. the external force Foxt applied on monomer 2 {£p/L = 1667) before the bifurcation for 
= Q.mip/L (□) and after = QAAip/L (o). c) Efficiency, 77, vs. Fcxt (same parameters as b). The solid line is a parabolic 

fit. 



to be balanced by the external torque N . The transition 
occurs when the bending radius reaches the length of the 
filament, R ^ L. This yields the critical torque 



TVc ~ -^sin^i 



(3) 



which is shown in Fig. ^ as a solid line and agrees very 
well with the numerical data. Neither hydrodynamical 
parameters nor temperature appear in Eq. (PI; the bi- 
furcation is thus purely elastic in origin and indeed can 
be obtained within free-draining or slender-body approx- 
imation (as will be published separately). Using the ro- 
tational mobility /.tjot ~ l/('7^^)j the critical angular 
frequency reads lOc ^ /^rot-^c ^ ^/iv^'^)- This thresh- 
old turns out to be much lower than for the continuous 
twirling-whirling transition of a rotating rod with tor- 
sional elasticity, w*"" ~ ^/(r/a^L^), which was obtained 
using linear analysis [l^. 

We define the translational and rotational mobilities of 
the filament as 1] 
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where Vz is the translational velocity downwards, fext 
the corresponding external force applied at the propeller 
base. In general, the mobilities depend on torque N 
and force fext in a nonlinear fashion. For perfectly stiff 
propellers, the mobility matrix is constant and symmet- 
rical, i.e. /irt = A^tr Q- However, in our case, due to flex- 
ibility, the propeller shape changes with external efforts 
and the matrix is asymmetric. The propulsion velocity 
along the rotation axis is plotted in Fig. as a func- 
tion of NL I £p for different persistence lengths for the z- 
moving case (i^cxt =0). At the transition, a jump in the 
propulsion velocity is observed and Vz is almost linear for 
N > Nc, revealing that the shape remains almost fixed 
in this range of torque values. Inset of Fig. |2K shows 
the ratio of Vz after and before the transition vs. the 
same ratio for lo. The variation is roughly linear, mean- 
ing that the jump of Vz is directly due to the jump in lo. 



To test our propeller under load, we applied an external 
force, foxt, which we define to be positive when it pushes 
against its natural swimming direction. Figure l^b shows 
the variation of Vz with Fcxt at two different torques just 
below and above the bifurcation. The laws are almost 
linear in both cases, meaning that fitt is almost indepen- 
dent of Fcxt ■ The efficiency of the power converter can be 
defined as the ratio of the propulsive power output and 
the rotary power input. 
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By inserting Eq. Q in Eq. ©, we obtain ?7(i^oxt) = 
-ititrNFcxt + Aitti^oxt)/(MrrA^^ + fl.tNFcxt)- We have 
checked that /irt is negligibly small. Hence, the efficiency 
becomes parabolic as a function of the external force as 
shown in Fig. The highest efficiency is obtained for 
-F'cxt = -fitiN/{2ntt) = -Fstan/2 and is only of the order 
of 1% after the transition and 3 times smaller before. 

Up to now we considered an isolated rotating filament 
to which external torque and forces were applied. In re- 
ality, the filament is attached to a base, and rotation of 
the filament is caused by some relative torque generated 
between base and filament. As for any self-propelling 
object in viscous solvent, the total force and torque add 
up to zero If the filament rotates counterclockwise 
at frequency w, then the base rotates clockwise with an 
angular velocity Q. Within our Stokes simulation, we 
model the base by three elastic arms of length La and 
persistence length which are coplanar and connected 
to the filament at one point (see Fig.|3K)- Hydrodynamic 
screening and lubrication are accurately accounted for on 
scales larger than the monomer radius a. The torque act- 
ing on the filament is balanced by a counter torque on 
the base. The resulting complex helical trajectories of 
this nanomachine make calculation of the mean propul- 
sion velocity difficult due to thermally caused reorienta- 
tion of the complete device. In the following, we therefore 
turn off the thermal noise (which plays a minor role for 
stiff polymers). Figure |21 shows the propulsion velocity 
Vp and the ratio of rotational speeds of the base and 
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FIG. 3: a) Sketch of the self-propelhng nanomachine whose 
base is made of 3 arms of length La- b) Propulsion velocity Vp 
(o) and ratio of the base, Q, and filament, u, angular velocities 
(■) against L/La for L = 30a, £p/L = 33, i^/ip = 10 and 
A'^ = 2000. The dashed line is the asymptotic velocity reached 
by an isolated filament. 

the filament, —fl/uj vs. the ratio L/La at fixed relative 
torque N = 2000 {L = 30a, ip/L = 33, i^/ip = 10). For 
base sizes comparable or larger than the filament length, 
L/La < 1 , the propulsion velocity is highly reduced com- 
pared to the case of an isolated filament (dashed line), 
mostly due to hydrodynamic friction caused by the base. 
But for decreasing base size, the propulsion velocity ap- 
proaches quite closely the velocity of an isolated filament, 
calculated according to the preceding model. Indeed, 
when it is very small, the counter-rotating base causes 
negligible hydrodynamic effects on the filament. At the 
same time, the base rotational speed fl (rescaled by w, 
see Fig. O increases when the base size decreases while 
remaining much smaller than the filament one. This is 
because, for such a torque much larger than the critical 
torque, the filament is close to the filament axis and thus 
gives little rotational resistance. It also means that only a 
small fraction of power input is invested into the base ro- 
tation. In conclusion, we show that a simple straight elas- 
tic filament can induce propulsion when rotated at one 



end. The interplay between elastic deformations and hy- 
drodynamic interactions results in a substantial directed 
thrust due to a subcritical dynamic bifurcation for con- 
trol parameters N / Nc > 1. The mobility /itr changes 
its sign when the torque is reversed, which implies a for- 
ward thrust whatever the sense of rotation. This propul- 
sion occurs even in the presence of an explicit base, being 
larger for small bases. This work provides a clue for the 
synthetic manufacture of biomimetic micropropeller, by 
using simple semiflexible polymers instead of rigid he- 
lices. This bifurcation could be experimentally investi- 
gated using a macroscopic scale model similar to the one 
developed by Powers group 0| . With their experimental 
setup, the bifurcation would occur for torques on the or- 
der of 0.05 Nm and angular velocities larger than 0.01 Hz, 
which are accessible. 

With regard to bacteria propulsion, flagella are heli- 
cal and their physics is therefore much more involved. 
However, flagellar motors generate sufficient torques for 
the nonlinear elastic phenomena discussed here: they are 
powered by a proton-motive force which yields torques 
10^ fcsT Jj5, 18J. Moreover, the flagellum length 

is i ~ 10 nm, leading to gp/L ~ 10^ lO'' and 

NL/£p ~ 0.1 1. Hence, the shape bifurcation dis- 
cussed in this model is probably biologically relevant 
and might be directly observed with straight biopoly- 
mers attached to flagellar motors. For a bacterial flagel- 
lum at physiological conditions, the diameter is roughly 
a ~ 20 nm and the stall forces are roughly 10 pN at the 
bifurcation. Critical angular and propulsion velocities 
(with no external force) then follow from our results as 

^ 2x10^ 10^ rad/s and T4c ^ 6 40 mm/s. These 

values are large but if we actually consider a thicker bun- 
dle formed by approximately 7 flagella |l3j |. we come up 
with smaller velocities which are comparable to the ones 
observed for E. coli {lo ~ 10** rad/s and 14 ~ 30 /im/s). 
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